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ABSTRACT 

Using Liverpool Telescope+RISE photometry we identify the 2.78 hour period binary star CSS 
Aim as a detached eclipsing double white dwarf binary with a 21,100K primary star and a 10,500K 
secondary star. This makes CSS 41177 only the second known eclipsing double white dwarf binary 
after NLTT 11748. The 2 minute long primary eclipse is 40% deep and the secondary eclipse 10% 
deep. From Gemini+GMOS spectroscopy we measure the radial velocities of both components of 
the binary from the Ha absorption line cores. These measurements, combined with the light curve 
information, yield white dwarf masses of Mi = 0.283 ± 0.064M© and M 2 = 0.274 ± 0.034M Q , making 
them both helium core white dwarfs. As an eclipsing, double-lined spectroscopic binary CSS 41177 is 
ideally suited to measuring precise, model-independent masses and radii. The two white dwarfs will 
merge in roughly 1.1 Gyr to form a single sdB star. 

Subject headings: binaries: eclipsing — stars: individual (CSS 41177) — white dwarfs 



1. INTRODUCTION 

There are of order 100-300 million close pairs of white 
dwarfs in the Galaxy. They are thought to be the domi- 
nant sour ce of backgroun d gravitational waves detectable 
by LISA (jHils et al.lll990t) . The loss of orbital angular 
momentum via gravitational radiation in these systems 
will eventually bring the two white dwarfs into contact 
with one another. Those that achieve this within a Hub- 
ble time are possible progenitors of A M CVn binaries, 
R Cr B stars and Type la supernovae (|Iben fc Tutukovl 
1984). The coalescence of two low-mass, helium core 
white dwarfs is also a possibl e formation cha nnel for the 
creation of sing le sdB stars ( Wcbbinkl Il98l : indeed it 
may be one of the mo re important sdB formation chan- 
nels (|Han et all 120031) . There are now 43 detached dou- 
ble w hite dwarf binaries known (Kilic et al. 2 0111: iMarshl 
120111 and references therein) but few of these yield the 
precise parameters needed to determine their eventual 
fate. 

The subject of this paper, CSS 41177 (SDSS 
J100559. 10+224932.2), was listed as a whit e dwarf plus 
main sequence binary by iDrake et al.l (|2010f) . They used 
a marginal 2MASS detection to infer the presence of 
a late-type companion though they note that a small 
faint object could also produce the observed transit sig- 
nal. We observed CSS 41177 as part of a monitoring 
campaign for period variations in white dwarf plus main 
sequence binaries. With a better sampled eclipse and 
detection of the secondary eclipse we realised that it 
was in fact a detached eclipsing double white dwarf bi- 
nary, making it only t he second known after NLTT 11748 
(jSteinfadt ct al. 20l3), but in this case the secondary star 
contributes 22% of the overall flux (compared to 3% in 



NLTT 11748). Follow up spectroscopy revealed that CSS 
41177 was also a double-lined spectroscopic binary. 

Here we present Liverpool Telescope (LT)+RISE pho- 
tometry and Gemini+GMOS spectroscopy of CSS 41177, 
and use these observations, and mass-radius relations, to 
determine the system parameters. We also discuss the 
evolution and eventual fate of the binary. 

2. OBSERVATIONS AND THEIR REDUCTION 
2.1. LT+RISE photometry 

Three primary eclipses and two secondary eclipses 
of CSS 41177 were observed using the LT and RISE 
([Steele et. al.1 120041 ) between February and April 2011. 
RISE is a high-speed frame tr ansfer CCD c amera with 
a single wideband V+R filter iSteele et al.l ([2008). All 
observations were made in 2 x 2 binning mode, giving a 
scale of 1.1 arcsec per pixel, and with exposure times of 
12-13 seconds. 

The raw data are automatically run through a pipeline 
that debiases, removes a scaled dark frame and flat-fields 
the data. The source flux was determined with aperture 
photometry using a variable aperture, whereby the radius 
of the aperture is scaled ac cording to the FWH M, using 
the ULTRACAM pipeline (jDhillon et al.ll2007D . Varia- 
tions in observing conditions were accounted for by de- 
termining the flux relative to the nearby star G 54-11 
(V=14.6). 

2.2. Gemini+GMOS spectroscopy 

CSS 41177 was observed using GMOS on the 8m Gem- 
ini Observatory North telescope ([Hook et al.l l2004) on 15 
March 2011. We used the R831 grating with a 0.5 arc- 
sec slit and 4x2 binning (spatial by spectral) to give a 
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Fig. 1. — Left: LT+RISE photometry of CSS 41177 folded on its 2.78 hour orbital period and binned by a factor of two. Right: 
unbinncd primary and secondary eclipse data. The secondary eclipse is shifted back by half a cycle and vertically by 0.3. 
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resolution of R ~ 4400 centered at Ha. We used an ex- 
posure time of 10 minutes, recording 17 spectra covering 
an entire orbit. 

The average bias and flat-field correction was car- 
ried out using the FIGARO package from STARLINK and 
nightly bias and tungs ten frames. We optimally ex- 
tracted each spectrum (Marsh 1989) and used a CuAr 
arc lamp to establish an accurate wavelength scale. The 
arc spectrum was fitted with a fifth-order polynomial to 
15 lines which gave an rms of 0.05 pixel. The individual 
spectra were normalised to the continuum level using a 
polynomial fit to the continuum regions and the spectra 
were placed onto a heliocentric wavelength scale. 

3. RESULTS 

3.1. Light curve model fitting 

The left hand panel of Figure ED shows the LT+RISE 
photometry of CSS 41177 folded on its 2.78 hour orbital 
period and binned by a factor of two. Since both eclipses 
are detected, the light curve information alone constrains 
the system inclination and the radii scaled by the binary 
separation. 

To measure the system parameters we fitted the light 
curve using a code written to produce models for the 
general case of binaries containing a white dwarf (see 
ICopperwheat et aLll2010l for details). The program sub- 
divides each star into small elements with a geometry 
fixed by its radius as measured along the direction of 
centres towards the other star. From an initial set of pa- 
rameters defined by the user, the code produces model 
light curves which are initially fitted to the RISE data 
using Markov chain Monte Carlo (MCMC) minimisation 
to produce a set of covariances. We then ran another 
MCMC minimisation using these covariances to define 
the parameter jumps. To produce the final parameters 
and their errors we followed t he procedures described in 
iCollier Cameron et al.l ()2007f ). The parameters needed 
to define the model were: the mass ratio, q = M 2 /Mi, 
the inclination, i, the radii scaled by the binary separa- 
tion, R\/a and TLij a, the surface brightness ratio, S1/S2, 
quadratic limb darkening coefficients for the both white 
dwarfs, the time of mid eclipse, To and the period, P. 

We corrected all the times to Bary centric Dynami- 
cal Time (TDB). We initially fitted each light curve 
individually in order to measure mid-eclipse times for 
the three primary eclipses recorded. The second eclipse 



TABLE 1 
Mid-eclipse times. 



Cycle 


Eclipse time 


U ncert 


number 


MJD(BTDB) 


MJD(BTDB) 


-171 


55599.087788 


0.000013 





55618.926463 


0.000016 


302 


55653.963116 


0.000020 



recorded was set to cycle number zero. Our measured 
eclipse times are shown in Table [TJ The zero-point of 
the ephemeris and the period were allowed to vary in 
our final fit. We kept the mass ratio fixed as 1.0 (the 
light curves of this well-detached system are indepen- 
dent of the mass ratio). We determined the quadratic 
limb darkening coefficients from a pair of white dwarf 
atmosphere models with temperatures of 21,100K and 
10,500K and log g = 7. 3 (see Section 15^1) fol ded through 
the RISE filter profile (jGansicke et al.lll995D . We deter- 
mined coefficients of a — 0.105 and b = 0.228 for the 
primary and a = 0.176 and b = 0.288 for the secondary 
for J(/i)//(l) = 1 — a(l — fi) — 6(1 — fi) 2 , where n is the 
cosine of the angle between the line of sight and the sur- 
face normal; we kept these values fixed. The inclination, 
scaled radii and surface brightness ratio were allowed to 
vary. 

The right hand panel of Figure Q] shows the best fit 
to the RISE photometry around the eclipses. The two 
eclipses constrain the system inclination to 89.2° ± 0.3°. 
The constraints on the two scaled radii are shown in 
the left hand panel of Figure [31 We also update the 
ephemeris to 

MJD(BTDB) = 55618.926447(9)+ 0.116 015 49(5)£, 
which is consistent with the ephemeris of IDrake et al.l 

(poTol) . 

3.2. Temperatures 

We used data from the Sloan Digital Sky Survey 
(SDSS) to determine the temperatures of the two white 
dwarfs in CSS 41177. Initially, we fitted its SDSS 
spectrum as a singl e white dwarf using the models of 
IKoester et all (|2005l ), which results in a T eSA = 21, 535 + 
214K and logg = 7.36 + 0.44. To determine the tempera- 
ture of both stars better, we fitted the SDSS ugriz mag- 
nitudes with a composite of two DA white dwarf mod- 
els, with the additional constraint of the surface bright- 



The Eclipsing white dwarf binary CSS 41177 




Fig. 2. — Left: Trailed spectrum of the H« line, the scale runs from black (50% of the continuum level) to white (90% of the continuum 
level). The absorption from the primary star is clearly visible. More subtly, absorption from the secondary star is visible moving in 
anti-phase with the primary stars absorption. Right: Same as left hand panel but with the primary stars contribution removed (black is 
-2% of the continuum, white is 4% of the continuum). The absorption from the secondary star is now visible. 



ness ratio Si/S 2 = 2.90 ± 0.23 in the LT/RISE band, 
as measured from the light curve. This fit results in 
T cfU = 21,100±800KandT cffi2 = 10, 500±500K, where 
the quoted errors are estimated, as they correlate with 
the uncertainties in the surface gravities. The best-fit 
composite model reproduces well the observed GALEX 
ultraviolet fluxes. 

3.3. Radial velocities 

The left hand panel of Figure [2] shows the trailed spec- 
trum of CSS 41177 around the Ha line. The non-LTE 
absorption core from the primary (hotter) star is clearly 
visible. There is also a weaker absorption component 
moving in anti-phase with the absorption from the pri- 
mary star originating from the secondary (cooler) star. 
No other features are visible in the GMOS spectra. 

Initially we fitted just the primary star's absorption. 
For each spectrum we fitted the line using a combination 
of a straight line plus a Gaussian. However, the radial 
velocity measured does not represent the true radial ve- 
locity of the primary star as the absorption from the 
secondary star causes a slight underestimation. Never- 
theless, we used these velocities to correct out the motion 
of the primary star. We then averaged the shifted spec- 
tra to obtain the rest-frame spectrum of the primary star 
which we then subtracted from all the observed spectra. 
The result is shown in the right hand panel of Figure [2j 
clearly showing the absorption from the secondary star. 

We fitted the absorption from the secondary star in 
the same way as the primary star. This too underes- 
timated the true velocity since the primary star's con- 
tribution had not been completely removed. However, 
both fits provide us with an initial starting point for K\ 
and K 2 , the radial velocity amplitudes of both stars. In 
order to measure accurate radial velocities both absorp- 
tion components need to be fitted simultaneously. We 
fit both lines together by simultaneously fitting all of the 
spectra. We used a combination of a straight line and 
Gaussians for each spectrum (including a broad Gaus- 
sian component to account for the wings of the primary 
star's absorption) and allowed the position of the Gaus- 
sians to change velocity according to 

V = 7 + if sin(27r0), 

for each star, where 7 is the offset in the line from its rest 
wavelength and <j) is the orbital phase of the spectrum. 
Due to the faintness of the secondary star's absorption we 



keep the offset 72 = 71 to reduce the number of degrees 
of freedom. All other parameters were initially set to 
those measured from fitting the lines individually. We 
fitted the data usi ng Levenberg-Marquardt minimisation 
(|Press et al.lll986T ). The resultant fit had a reduced x 2 of 
1.02. Wefindifi = 177±3kms" 1 , K 2 = 181±20kms~ 1 
and 7 = 109 ± 2 kms -1 . These values imply that CSS 
41177 is an equal mass binary (or very close to, given the 
uncertainty on K2). 

3.4. System Parameters 

The RISE photometry constrains i, R\/a and R2/a. 
K\ and K 2 are constrained from the spectroscopy, there- 
fore we can determine the masses and radii of both white 
dwarfs using Kepler's third law, 



G{M 1 + M a ) 



4tt2 

p2 ' 



(1) 



where P is the orbital period, given the orbital separation 
from 



2tt 

— a sin i = Ki + K 2 



(2) 



We determined the masses and radii for each model pro- 
duced in the MCMC analysis of the light curve using 
the measured radial velocities and Kepler's laws. The 
right hand panel of Figure [3] shows the regions of al- 
lowed masses and radii for both white dwarfs based 
only upon the photometric and spectroscopic constraints. 
Also shown are models of helium core white dwarfs 
with hydrogen layer fractional masses of Mh/M, = 
10" 4 (solid lines) and Mh/M* = 10" 8 (dashed lines) 
from iBenvenuto fc Althausl (|1998l ) for temperatures of 
10,500K (gray) and 21,100K (black) which, given the 
large uncertainties in our measurements, are consistent 
with the measured masses and radii. The final system 
parameters are listed in Table [2j 

4. CONCLUSIONS 

We have discovered that CSS 41177 is a detached, 
eclipsing, double white dwarf binary containing two low- 
mass helium core white dwarfs with masses of Mi = 
0.283 ±0.064M o and M 2 = 0.274 ±O.O34M and radii of 
Ri = 0.0210 ± O.OO26R and R 2 = 0.0174 ± 0.0031R©. 
The hotter white dwarf has a temperature of 21, 100 ± 
800K and the cooler white dwarf has a temperature of 
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Fig. 3. — Left: The contours show the 68 and 95 percentile regions of the scaled radii R\/a and Ri/a, where a is the orbital separation, 
constrained by the RISE photometry. Right: The contours show the 68 and 95 percentile regions of the masses and radii for the hot (black 
lines) and cool (gray lines) white dwarfs. Also shown are mass-radius relations for helium core white dwarfs with hydrogen layer fractional 
masses of M H /M, = 10~ 4 (solid lines) and M H /M» = 10 -8 (dashed lines) from IBenvenuto fc Althausl (119981) for temperatures of 10,500K 
and 21,100K. The confidence regions shown here are for two parameters jointly and should not be directly compared to the uncertainties 
in Table [2] which are more closely related to single parameter confidence intervals. 



TABLE 2 

Stellar and binary parameters for CSS 41177. 



Parameter 


Value 


Parameter 


Value 


RA 


10:05:59.11 


K L (kms" 1 ) 


177 ±3 


Dec 


+22:49:32.3 


K 2 (kms" 1 ) 


181 ± 20 


GALEX fuv 


16.72 ±0.04 


Teff.l (K) 


21, 100 ± 800 


GALEX nuv 


16.98 ±0.03 


T eff ,2 (K) 


10, 500 ± 500 


u 


17.317 ±0.010 


Mi (Mq) 


0.283 ± 0.064 


9 


17.266 ± 0.005 


M 2 (Mq) 


0.274 ± 0.034 


r 


17.612 ±0.006 


Ri (R©) 


0.0210 ± 0.0026 


i 


17.899 ± 0.008 


R 2 (R ) 


0.0174 ± 0.0031 




18.151 ±0.024 


a (R ) 


0.821 ± 0.048 


d (pc) 


350 ± 13 


i (deg) 


89.2 ±0.3 


Porb (days) 


0.1160156(1) 







10, 500±500K placing it on the red edge of the instability 
strip. 

CSS 41177 is ideally set up to measure precise model- 
independent masses and radii for both white dwarfs and 
thus can be used to test white dwarf mass-radius rela- 
tions. To date only one white dwarf has had its mass and 
radi us independently m easured to high precision (NN 
Ser. IParsons et~al]l2ninl) . therefore CSS 41177 offers two 
additional white dwarfs which can test mass-radius rela- 
tions at the lower mass range. 

Since the white dwarfs in CSS 41177 have roughly 
equal masses they will eventually undergo a merger re- 
sulting in the fo rmation of a single ~ O.6M0 sdB star 
(jHan et all [20031 ). This will occur in roughly 1.1 Gyr 
due to the loss of orbital angular momentum via gravi- 
tational radiation. The regular eclipse times from CSS 



41177 should make it an excellent timing standard for 
optical astronomy, provided there are no unseen com- 
panions. 



Balmcr/Lyman lines in the models were calcu- 
lated with the modified St ark broadening profiles of 
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